OBSERVATIONAL CONSEQUENCES OF BARYONIC GASEOUS DARK MATTER 



Yu. A. SHCHEKINOV 

Department of Physics, Rostov State University, Rostov on Don, Russia 

yus@rsussl . rnd. runnet. ru 



Abstract. Possible observational consequences of dark matter in the Galaxy in the 
form of dense molecular gas clouds - clumpuscules of masses M c ~ 10 -3 M Q and 
radii R c ~ 3 x 10 13 cm - are considered. Recent models of the extreme scattering 
events - refraction of radio-waves from quasars in dense plasma clumps in the Galactic 
halo - definitely show on such clouds as possible dark matter candidate. We argue that 
collisions of such clumpuscules are quite frequent: around 1—10 Mq a year can be ejected 
in the interstellar medium due to collisions. Optical continuum and 21 cm emissions 
from post-collisional gas are found to be observable. We show that clumpuscules can 
form around O stars HII regions of sizes R ~ 30 pc and emission measure EM ~ 20 
cm -6 pc, and can also be observable in H a emission. Evaporation of clumpuscules by 
external ionizing radiation can be a substantial mass source. From requirement that the 
total mass input on the Hubble time cannot exceed the luminous mass in the Galaxy 
typical radius of clouds is constrained as R c < 3.5 x 10 12 cm, and their contribution to 
surface mass density as 50 M Q pc -2 . It is argued that dissipation of kinetic energy of 
the gas ejected by such clouds can be an efficient energy source for the Galactic halo. 
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1. Introduction 



Baryon dark matter has been supported observationally when microlensing events by 
massive compact halo objects (MACHOs) were detected [1—3]. The nature of MACHOs remains 
still unclear, although from general point of view stars of small masses, M ~ 1O _3 M (brown 
dwarfs, called often "jupiters") can provide gravitational lensing with observed characteristics, 
and on the other hand, being below the limit of hydrogen burning, they have too low luminosity 
to be seen. It is quite possible that microlensing events with observed frequency can be connected 
also with red dwarfs (M £ 0.08 M Q ) with total mass exceeding the mass of luminous matter in 
the Galaxy [4]. 

Recently dense atomic and molecular clouds, confined by external pressure or self- 
gravity, are discussed as an alternative possible reservoir of dark matter [5-14]. Contrary to 
brown dwarfs these clouds are assumingly cold (from 3 K [6] to 10 K [14]) and have negligible 
proper luminosity and are unseen even in infrared, their column densities are predicted to reach 
~ 3 x 10 23 [11, 12] to 10 25 cm -2 [13], and they are optically thick in 21 cm. From this point of 
view dark matter gas clouds are notoriously hard to detect. 

In [7] 7-radiation from interaction of cosmic rays with nuclei of dark matter molecular 
clouds has been estimated and suggested to test the presence of such clouds in the galactic halo. 
Very recently statistical evidences for a large scale anisotropic excess in 7-rays were presented 
[15], with the flux above 1 GeV close to the one predicted in [7]. This can be considered in 
support of dark mass in the form of dense cold gas. 

Additional evidences for the existence of dark matter gas clouds are connected with 
optical (non-gravitational) lensing of radio emission from quasars passing through atmospheres 
of ionized gas surrounding such clouds. Ten years ago dramatic short variations of fluxes (of 
several weeks to months duration) from compact radio quasars were discovered at wavelengths 
3.7 and 11 cm - the so-called extreme scattering events [16]. The origin of these variations 
was attributed to refraction by dense plasma clouds moving perpendicularly to the line of sight 
[16, 17]. From variations of radio fluxes lensing clouds are inferred to have characteristic sizes 
of a few AU and electron densities of n e ~ 10 3 cm -3 . The internal gas pressure n e T e £ 10 7 
K cm -3 exceeds pressure in the surrounding gas by 3-4 orders of magnitude, and therefore 
refracting plasma clouds must be transient structures of several months life time, which is far 
from being enough to explain observed frequency of the extreme scattering events [16, 17]. Only 
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recently this difficulty has found a solution [18] with the assumption that lensing plasma can be 
connected with extended HII atmospheres around cool dense self-gravitating molecular clouds 
as was suggested in [9] . In the framework of this model high interior pressure provides steady 
outflow of ionized gas from the HII cloud atmosphere, and thus the density distribution in 
lensing plasma regions is kept time-independent. The individual cloud masses are high enough 
to support photoionization mass outflow on galactic time scales. In [18] these clouds were 
identified with the dense molecular clumpuscules described in [5, 6]. This allowed the authors 
to reproduce with good accuracy the dual- wavelength (3.7 and 11 cm) light-curve of the quasar 
0954+658, while all the previous models had difficulties with reproducing this. Quite recently, 
the idea that such dense molecular clouds contribute significantly to dark matter has obtained 
further support: it was shown that they act as gaseous (non-gravitational) lenses for stellar 
light, which in turn can result in the observed excess of microlensing events, compared to the 
expected from stars and stellar remnants toward the LMC [19]. In this paper we consider 
possible dynamical and observational consequences of the existence of dark matter in the form 
of dense molecular clouds. (Through the paper we will refer these clouds to as DM-clouds, or 
following [13], clumpuscules). 

In §2 we will discuss general properties of the population of DM-clouds in the Galaxy 
which can be inferred from the results of [18]; in §3 we describe qualitatively dynamics of 
collisions of DM-clouds; in §4 we estimate optical emission from colliding clumpuscules; in §5 we 
estimate 21 cm emission from clumpuscules and their debris after collisions; in §6 we consider 
H a emission from ionized atmospheres of DM-clouds; in §7 possible implications for the Galactic 
halo gas is discussed; §8 contains discussions of the results and conclusions. 

2. General properties of molecular DM-clouds. 

In [5, 6] molecular DM-clouds were assumed to have high enough total mass (about one 
order of magnitude in excess of the mass of visible matter), quite low proper luminosity, and life 
time comparable with the Hubble time. It follows then that individual masses of these objects 
must be close to the Jovian mass, M c ~ 10~ 3 M , and densities, temperatures and sizes such 
to make them marginally gravitationally stable: n c ~ 10 11 cm -3 , T c ~ 3 K, R c ~ 10 14 cm. The 
interpretation of the extreme scattering events suggested in [18] leads to very close value for the 
radius of a DM-cloud: the duration of these events around 2 months suggests that R c ^ 3 x 10 14 
cm for the transverse velocity of DM-clouds restricted from above by the escape velocity for the 
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Galaxy 500 km s _1 . The covering factor of the clouds is estimated from the flux monitoring 
data [16] as / ~ 5 x 10~ 3 . It allows, in turn, to estimate the contribution of DM-clouds to the 
total surface density of the Galaxy as [18]: £ ~ fM c /irRl, where M c is found from hydrostatic 
equilibrium M c ~ kT c R c /Gm p . It gives for T c > 3 K the surface density £ > 10 2 M Q pc~ 2 [18], 
which is close to the value estimated from dynamical arguments [20]. 

The characteristic free-path length of DM-clouds with respect to cloud-cloud collisions 
is of r //, where r is the characteristic size of the region occupied by clouds. For r ~ 10 kpc, 
one gets for characteristic time between collisions ^ 10 10 yr. Cloud destruction by the external 
UV radiation is slow enough with similar characteristic time, due to their high masses and 
densities. Moreover, contact interactions of DM-clouds with interstellar clouds are inefficient 
because of high differences in densities, and thus one can conclude (see [18]) that the population 
of DM-clouds does not change significantly over the Hubble time. At the same time, even so 
inefficient destruction of DM-clouds can result in dynamically important consequences because 
their total mass in the Galaxy is large. 

We will assume for specifity the distribution of DM-clouds in the form 



which fits the singular isothermal sphere at r > vq (through the paper rg is regarded as a free 
parameter); here Ao is the number of clouds per unit volume in the galactic center, r is the 
galactocentric radius. The collision rate per unit volume is then 



where all clouds are assumed to have equal characteristic velocity v c independent of their position 
in the Galaxy. In what follows we accept in numerical estimates R c = 3 x 10 13 cm and v c = 250 
km s -1 - two times less than the escape velocity - unless other values specified. The total rate 
of collisions in the Galaxy is 
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which is connected with the characteristic covering factor of DM-clouds / ~ tiNqR^tq by equa- 
tion 



4 r ° Vc f2 (a\ 

v G ~ -jgf ■ (4) 

For adopted values of R c , v c and / this gives vq ~ 10 _26 r s _1 , or vq ~ (10 3 — 10 4 ) yr _1 
for r = 1 kpc and r = 10 kpc, respectively; the total number of clumpuscules in the Galaxy 
can reach N ~ 10 13 — 10 15 . (Note, that for given characteristic length of clouds distribution 
r and fixed total number of clouds in the Galaxy N, the surface density from clumpuscules is 
proportional to their radius £ ~ kT c NR c /-KGm p rQ, however, for fixed current covering factor 
it is inversely proportional to R c : £ oc fRj 1 .) Therefore, the amount of material of molecular 
clumpuscules processed by collisional shocks is ~ vqM c ~ 1 — 10 M Q yr _1 . This value is 
comparable or exceeds typical rates of matter processing in such galactic events as star formation, 
heating of the interstellar gas by supernovae shocks, accretion of galactic halo gas to the plane. 
One can expect thus that collisions of DM-clouds have important consequences for the interstellar 
gas and might be observed. 

3. Collisional dynamics of molecular DM-clouds. 

In this Section we describe qualitatively dynamics of collisions of DM-clouds. Detailed 
description of numerical results will be given elsewhere. Let us consider a head-on collision of two 
equal clouds, and neglect self-gravity effects. In the center-of-mass frame a contact discontinuity 
forms in the symmetry plane, from which two shocks propagate outward. Characteristic duration 
of the collision, i.e. time needed for the shock to reach outside boundary of the cloud, is r s ~ 
R c / V c ~ 0.4 yr for adopted parameters. Due to high density behind the shock, n s = 4n c ~ 
4 x 10 12 cm -3 , radiation cooling time in bremsstrahlung processes is only tr ~ 3 s. However, 
hard radiation emitted by shocked gas is absorbed on scales ~ 10 6 cm, and after time interval 
~ 10 5 s, when the shock wave covers around 1 - 2 % of the cloud mass, it becomes ionized 
x ~ 1 and Thomson optical depth increases to tt ~ 20. This circumstance allows us to neglect 
radiative cooling. 

Due to high postshock pressure, P s ~ p c v 2 c ~ 10 3 erg cm -3 (~ 10 19 K cm -3 ), a strong 
axial gas outflow along the symmetry plane forms immediately after collision, which decreases 
pressure and temperature behind the shock, and moreover it scatters a substantial mass fraction 
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of clouds during the collision t ~ t s . In [21] this effect was described first and then confirmed 
in a sequence of papers [22-26]. The fraction of scattered mass for diffuse interstellar clouds 
dominated by radiative cooling was estimated in [21] as 25 %. In radiationless case this fraction 
can be somewhat larger due to higher pressure. (One should stress, however, an approximate 
character of these estimates - recent simulations with high precision [26] have demonstrated that 
at late stages after collision, t > 10 — 15 dynamical times, Rayleigh- Taylor instability develops to 
form multiple clumps and filaments, and violates strongly the integrity of post-collisional clouds. 
Thus, the fraction of escaped gas can be higher.) 

After the shock wave reaches external boundary of the cloud, a rarefaction wave starts 
to propagate inward resulting in expansion of the cloud with velocity equal to the local sound 
speed (see [21]). At t ^ t s the rarefaction wave reaches the symmetry plane and reflects. When 
the reflected wave reaches external boundary of the cloud, gas pressure in outer regions relaxes 
already to external pressure, and thus the secondary reflected wave is a compression wave. These 
wave motions result in separation of clouds with final pressure relaxed to the external value. 
Since viscosity behind shock waves works to increase entropy, cloud collisions are dissipative even 
in radiationless case. This means that after separation of clouds pressure equilibrium establishes 
at gas density less than the initial one: the new equilibrium corresponds to the Poisson adiabat 
with entropy larger than the initial value. 

Self-gravity effects on collisional dynamics are weak, and are only important for unper- 
turbed clouds (i.e. before collision), as they are assumed to be marginally stable [5, 6]. After 
collision thermal energy of clouds becomes larger than their gravitational energy due to dis- 
sipative increase of entropy at the shock, and separated clouds expand freely up to pressure 
equilibrium with ambient medium. At late stages clouds can get transparent to external heating 
sources, and their subsequent expansion settles in an isothermal regime, and finally they mix 
with surrounding gas. 

An important qualitative difference of shocked gas in colliding clumpuscules from that 
described in [21 - 26] is connected with the fact that due to high relative velocity of clumpuscules 
the postshock gas is dominated by radiation pressure. As mentioned above, Thomson optical 
depth of the shocked gas reaches unity quickly, and thus gas and radiation behind the shock set 
in equilibrium. The postshock pressure is then 
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P=KpT+-a SB T\ (4) 

fj/ O 

where 1Z is the gas constant, a SB [erg cm -3 K -4 ], the Stefan-Boltzmann constant. On the other 
hand, P = pqD 2 (1 — rj), where r] -1 = p/po is the compression factor, D, the shock velocity 
in the center-of-mass frame. It can be readily shown that the solution of eq. (4) at v c = 250 
km s _1 is T ~ 2.8 x 10 4 K, and the contribution of gas pressure [first term in eq. (4)] is 
~ 0.01 for n c = 10 12 cm -3 . Fractional ionization (determined by Saha equation) is x = 1 (more 
precisely 1 — x ~ 6 x 10 -7 ). Due to radiation diffusion the shock front is smoothed (see [27]) over 
scales of several free-path lengths of photons, in our case ~ 10 _2 i? c . From this point of view, 
qualitative dynamics of shocks and subsequent separations of clouds after collision is similar to 
that described in [21]. The sound speed behind the shock is 

s 3 pV2 ' 
and for adopted parameters c s ~ 0.7 x 10 7 cm s _1 . 

4. Optical emission. 

The postshock temperature is 

for v c ^ 15 km s _1 (for n c ~ 10 12 cm -3 ), at lower velocities the contribution of radiation to 
pressure behind shock is smaller than gas pressure, and temperature is determined by 

T s ~2.5xl0 6 ( ^ ^ 2 K . 

V250 kms- 1 / 

Thomson optical depth after ionization of the preshock gas reaches tt ~ aTR c n c ~ 20. The 
luminosity of the colliding DM-clouds is determined by diffusion of photons from the hot interior 
and is equal to 

L h ~ a SB T^cT 4 S, 

where S ~ 2-7ri? 2 is the emitting area. For adopted parameters it gives ~ 10 38 erg s _1 . 
The duration of hot stages is determined by gas expansion after shock waves reach external 
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boundaries of the clouds and the rarefaction waves start to propagate inward: td ~ R c /3c s ~ 
2 x 10 6 s. Thus, the total energy loss by radiation is Lhtd ~ 2 x 10 44 erg, i.e. about 30 % of the 
total kinetic energy of clouds. 

Observationally colliding DM-clouds at these stages are short-lived transients (with 
the life time t ~ t d ~ 10 d ), with temperature and luminosity typical for massive stars. The 
total number of such objects in the Galaxy is vcta ~ 80(r /l kpc). Only small their fraction 
~ 27r _1 arctg(r/r ), is located inside a region with galactocentric radius less than r, and can be 
obscured by dust in the galactic central parts: this fraction inside r < 1 kpc is 0.5 for rp = 1 
kpc {i.e. the number of unobscured objects is ~ 40), and only 0.06 for r$ = 10 kpc. Their 
characteristics (for v c > 100 km s _1 and n c ~ 10 12 cm -3 ) scales as the total number of clouds 
at hot stage 

Nh = — (——^\ X ^ / 2 "77" — 80 ( r ° \ ( Rc 



vr V v J R c V 1 k P c / V 3 - 1013 cm 

cloud luminosity 
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Lh = T J ^^~7.5xl0 37 -— or;n1 — ergs" 

1 — rj cjt \ 3 • 10 id cm y V 250 km s 1 / 



temperature 

1/4 



V 2 PqVc 
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2.8 x 10 4 ( -—^ ? I ( I K, 

V10 12 cm- 3 / V250 kms- 1 / 



here rj ~ 1/6. Possible detection of such objects (or contrary, non-detection at given threshold) 
might allow to infer (or restrict) characteristics of DM-clouds and estimate their contribution 
to the surface density of the Galaxy. 

5. Emission in 21 cm. 

Estimates of emission in 21 cm after collision of clouds can be obtained assuming each 
cloud to expand isotropically after the primary shock reaches the external surface of the cloud. 
At such an assumption, initial stages can be described as an expansion with 7 = 4/3 

\ 4 /3 



P = Pl 

.Pi. 

where subscript 1 refers to variables at the moment of maximal compression: Pi = 4ct sb T 4 /3, 
Ti = 2.8 x 10 4 K, pi = p c /r). At stages when Thomson optical depth is large, tt > 1, both 
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radiation and gas temperatures vary as T oc i? . For a quasi-spherical expansion with gas 
density varying as p oc i? -3 , both gas (P g ) and radiation (P 7 ) pressure vary as oc i? -4 , thus the 
initial ratio Pg/P-y — 0.05 is kept up to stages when electron density decreases substantially and 
tt becomes less than one. It can be readily shown that tt = 1 is reached when temperature 
decreases to T 2 = 5.8 x 10 3 K. At this moment n 2 = 0.016 n\ ~ 10 11 cm -3 , x = 0.12, and column 
density of neutral hydrogen is N 2 {HI) ~ 6 x 10 23 cm -2 . Note, that hydrodynamical instabilities 
in post-collisional gas, found in recent numerical simulations, produce non-homogeneous density 
distribution in clouds and their fragments, however, the density contrast between different re- 
gions does not exceed factor of 3-4 at times t > 50 — 60 dynamical times [26] , which in our case 
corresponds to t ~ 10 6 — 10 7 s. This means, that our description based on the assumption of 
smooth structure of post-collisional clouds, gives correct estimates by the order of magnitude. 

Subsequent expansion is adiabatic with 7 = 5/3, and at 713 = 10 6 cm -3 gas temperature 
is T3 ~ 3 K. The duration of this stage is ~ R^/c s ~ 2 x 10 11 s, where R3 ~ (ni/n^) 1 / 3 /^ is 
the characteristic size of the cloud corresponding to density n = n^: for adopted numbers 
R 3 ~ 10 15 cm. Note, that possible existence of such cold cloudlets in the gaseous galactic 
halo is argued in [28] from analysis of turbulent motions in the neutral halo gas. Further the 
cloud expands isothermally with temperature T ~ 3 K, which is supported by deactivating 
inelastic collisions of hydrogen atoms with electrons H(F = 1) + e — ► H(F = 0) + e' , with a 
rate qio = 3 x 10 -11 T 1 / 3 cm 3 s _1 [29] - the corresponding characteristic time iio ~ 2 x 10 4 
s. Isothermal phase stops when gas density decreases to 114 ~ 3 x 10 3 cm -3 and the cloud 
comes to pressure equilibrium with the ambient interstellar gas. At this state heating by the 
background X-ray emission gets important with characteristic time tx ~ 4 x 10 11 s (the heating 
rate Tx — 10~ 27 erg s _1 H _1 is taken from [30]). Subsequently, the cloud expands slowly 
with charateristic time tx, and becomes transparent to 21 cm when HI column density falls to 
N(HI) ~ 10 19 cm -2 - it occurs at t = tui ~ 10 13 s when n = ~ 10 3 cm -3 . Apparently, 
at t > tni the cloud mixes with the interstellar gas. The total number of such clouds in the 
Galaxy is estimated as v G t H i ~ 3 x (10 8 - 10 9 ), and their total mass as 3 x (10 5 - 10 6 ) M Q 
for r = 1 — 10 kpc, respectively. A distinctive feature of these clouds is their high velocity 
dispersion which is close to the dispersion of DM-clouds v c (for DM-clouds with Maxwellian 
velocity distribution the mean center-of-mass velocity for pairs of clouds is ~ ^/3/2v c ) - much 
larger than the local velocity dispersion of the interstellar HI. This circumstance can be used for 
identification of such clouds. It should be mentioned, however, that when cloud fragments and 
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debris mix with surrounding gas, a fraction of their kinetic energy converges to thermal energy 
due to viscosity, and thus a fraction of their mass will be heated up to high temperatures and 
thus will be unseen in 21 cm. To evaluate this effect, note that drag force for a cloud moving 
through interstellar gas is proportional to ram pressure oc piV 2 Rl, even for subsonic clouds [31]. 
In these conditions the length scale for velocity decay is id ~ SR C , where 5 is the density contrast 
between cloud and intercloud gas [32]. At stages when fragments become transparent to 21 cm 
gas density is of order ~ 10 3 cm -3 and radius R4 ~ 10 16 cm. Thus, cloud velocity decays 
on scales Id ~ 30 — 300 pc for density of ambient gas rtj ~ 0.1 — 0.01 cm -3 , respectively. The 
corresponding characteristic time ta ~ 3 x 10 12-13 s is comparable to tui- Therefore, although 
dynamics of cloud fragments at late stages can depend on environment crucially, time interval 
for fragments to be seen in 21 cm is close to tui- At these stages, before mixing with the ambient 
gas, fragmets form a network with covering factor of ~ 0.01 — 0.3 and resemble the tiny-scale 
atomic structures described in [33]. 

DM-clouds of smaller radii, i? c ^3x 10 12 cm, have small Thomson optical depth and 
due to putchy density distribution generated after collision by hydro dynamical instabilities, [26], 
they are apparently transparent. In this pointed in [34], after collison they quickly lose 

their energy radiatively and escape the hot phase. However, subsequent expansion, starting 
from stages when temperature has fallen to several thousands, is qualitatively similar to the 
described above, though it is less rapid since at given pressure clouds have larger densities. 

6. H a emission from isolated DM-clouds. 

By isolated clouds we mean clouds escaping collisions. Gas temperature in extended 
ionized gas around DM-clouds is of order T ~ 10 4 K, which for electron density n e ~ 10 3 cm -3 
produces pressure by 3-4 orders of magnitude in excess of typical interstellar pressure. In such 
conditions ionized gas expands freely with the velocity no = ^2/7 — 100, where cq ~ 10 km s _1 
is the sound speed. Therefore, ionized atmospheres of DM-clouds are in steady outflow regime 
with density and velocity distributions described by the equations (see, e.g. [35]) 

pur 2 = PqUqR 2 c , 

^ + ^\ n A = ^ 
2 p p 2 ' 
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which give for density distribution the solution 
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The equation of ionization equilibrium 
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with n e oc p taken from equation (6) determines the boundary value n e0 at r = R c as 



n e o - 



2J 



where ( = r/R c , and J [cm -2 s _1 ] is the intensity of ionizing photons at Lyman edge, a r , the 
recombination rate. The corresponding emission measure is 



EM = J nidr = 0.236 n z e0 R c = 1.8 cm" 6 pc, 

Rc 

here we adopted J obtained in [36] at high galactic latitude J = Jo ~ 10 6 cm -2 s _1 sr _1 , or 
the intensity of H a photons of ~ 1 R (= 10 6 /47r phot cm -2 s _1 sr _1 ). This number is by an 
order of magnitude higher than obtained in [37] for high-velocity HI clouds. However, because 
of small angular sizes H a fluxes from distant DM-clouds at d > lpc with angles < 0'.2 are 
too weak to be detected at present. At the same time, clouds at distances d ~ 0.03 pc with 
angular sizes A(p ~ 7' could be seen in H a . For adopted parameters and vq = 1 kpc, a volume 
d 3 ~ (0.03) 3 pc 3 in solar vicinity contains in average one DM-cloud; for rg = 10 one cloud is 
contained in a volume with radius 0.07 pc - at such distance cloud angular size is only ~ 3', 
and H a flux decreases by factor of 5. 

Much larger emission measures have clumpuscules located in close vicinity of hot stars 
where the intensity of ionizing photons is higher than the background value J » Jo- For 07 
stars this condition holds at r < 30 pc. The emission measure increases proportionally to J/ Jo, 
and thus a DM-cloud at distance ~ 1 pc from hot star can be seen as an ultracompact H a 
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nebula as pointed in [12, 18]. However, because of small angular sizes (~ 0".l for a cloud at 
~ 100 pc from Sun), expected fluxes are too low: Fn a ~ 3 x 10~ 3 phot cm -2 s _1 from a single 
cloud within 1 pc near an 07 star. In general, the total number of H a photons emitted by a 
single cloud is 

CR 2 

L(r) ~ 0.025^ phot s _1 , (7) 

where C [phot s _1 ] is the number of ionizing photons from a star, r, the distance to the star. A 
distance from the star where the background UV flux is equal to the flux produced by star is 



R r , 



4vrJ ' 

R m ~ 30 pc for an 07 star. Total number of H a photons emitted by all DM-clouds inside 
r < R m is then determined as L m = 4nAf J^ m L(r)r 2 dr and is equal 

L m = 0.3^=^ phot s- 1 . 
\/47rJ 

For adopted numerical values this gives ~ 7 x 10 44 phot s _1 , however the corresponding surface 
brightness is extremely small - the equivalent emission measure is only EM ~ 0.0025 cm -6 pc. 

One of possible observational manifestations of DM-clouds is connected with the fact 
that outflowing ionized gas fills the space around an O star, and if the outflow rate is high 
enough it can support high density of electrons and can be observed as a luminous halo. In 
order to estimate this effect let us write the approximate equation of mass balance at distance 
r from an ionizing star as 

^(pu) =ATM c (r)r 2 , (8) 

where M c = In R 2 m p n e Q(r)uQ is the mass loss rate by a single cloud (factor 2ir accounts that 
only one side of the cloud is ionized by the star), n e0 (r) = 4.8 x J 1 / 2 ^ 1 ^ 2 , (thus n e0 oc r _1 in 
virtue of J oc r~ 2 ). Note, that a volume with r ~ 10 pc contains ~ 10 6 to ~ 10 7 clouds, and this 
substantiate the assumption of steady state mass and energy balance. In this case, integration 
over r gives 

m p nu = NM c (r)r, 
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or finally 



n = 8.5 x 10 6 C^R^N^. (9) 

u 

To complete description, Bernoulli and energy equations must be written. As we will 
show below, gas pressure in regions filled by gas from HH-atmospheres is much higher than 
pressure in surrounding interstellar gas, and thus one can assume the outflow velocity u to be 
independent of r and equal with the factor 1/2/7 — 1 to the sound speed, which in turn is 
determined by energy equation. It is readily seen that in this case n =const. Energy budget of 
the gas is determined by heating radiation from central star and by dissipation of high-velocity 
turbulent motion of the gas lost by DM-clouds from one side, and radiation cooling from the 
other. The mean velocity of the gas lost by clouds relative to intercloud gas is approximately 
equal to clouds velocity v c , and thus the heat released due to dissipation of kinetic energy of 
evaporated gas can be considerable, ~ 10 14 erg g _1 (the corresponding temperature T ~ 10 7 
K). At such assumptions energy equation can be written as 

H uv (T)n 2 + \MM c vl = A(T)n 2 , (10) 

where Huy(T)n 2 is the heating rate by UV radiation from central star, Huy(T) = a r (T)eo, eo> 
the mean energy of photoelectrons (see [29, 38]), A(T), the radiation loss rate. Assuming DM- 
clouds to be of primordial chemical composition [6, 13], one can accept that radiation cooling in 
interval T = 10 4 — 10 6 K is determined by recombinations and equal to A(T) ~ 9.5 x 10 _17 T -3 / 2 
erg cm 3 s" 1 at T = 10 4 - 6 x 10 4 K and ~ 3 x lO" 21 ^ 1 / 2 at T = 6 x 10 4 - 10 6 K (see [29]). 
Eliminating n from (9) and (10) and assuming u = ^2/7 — lc s , we find that temperature varies 
from T ~ 10 4 K at r = 1 pc to T = 3 x 10 4 K at r = 30 pc. Therefore, one can accept in 
estimates T = 10 4 K in the whole range of r, which gives u ~ 20 km s _1 . Substituting u in (9) 
we arrive at 

n ~ 4.2 x 10 6 £ 1/2 i? 3/2 AA, (11) 

and 

EM ~ 1.2 x W- 6 CR 2 C M 4/3 cm" 6 pc, (12) 
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which gives n ~ 0.65 cm 3 and EM ~ 20 for adopted parameters and r = 10 kpc. The surface 
brightness in H a is large enough i# a ^ 8 R, and thus such HII halos can be observed. 

7. Contribution of DM-clouds to energy balance of the Galactic corona. 

DM-clouds spend most time in the halo, where their ionized atmospheres are supported 
by the background UV radiation [18]. The mass loss rate by a single cloud is at such conditions 



M c = 4Trm p J^u R 3 J 2 , (13) 
V QV 

or numerically M c ~ 10~ 13 M Q yr _1 . The evaporating gas loses its relative velocity due to 
viscous forces on characteristic length of Id ~ SR C (see §5), where radius of ionized atmosphere 
is accepted to be of order of cloud radius. For mean densities of halo gas ~ 10~ 3 cm -3 and of 
ionized atmosphere n e ~ 10 3 cm -3 this gives id ~ 10 6 R c ~ 10 pc which is much less than typical 
galactic scales. Thus, the energy input rate from kinetic energy of evaporating atmospheres is 
order of H c ~ M c Nv 2 c /2 ~ 2.5 x 10~ 25 erg cm -3 s -1 , which gives f H c dz ~ 2.5 x 10~ 3 erg 
cm -2 s _1 per unit area (here N ~ //7ri? 2 r , see §2, R c = 3 x 10 13 cm and r = 10 kpc are 
substituted). For comparison, the total energy input rate from SNe II is ~ 2 x 10~ 4 erg cm -2 
s _1 for SNe II galactic rate 1/30 yr _1 . The radiation loss rate for gas with temperature T ~ 10 6 
K and density n ~ 10~ 3 cm -3 is order of A(T)ra 2 ~ 3 x 10~ 29 erg cm -3 s _1 . In principle, 
observations show that the halo gas has extremely non-homogeneous density and temperature 
distributions. In this case, one can expect the total rate of radiation energy losses to be higher 
than the above estimate. In particular, in regions with n ~ 10 -2 cm -3 and T ~ 10 5 K radiation 
takes away ~ 10 -25 erg cm -3 s _1 . However, energy input rate from evaporating clumpuscules 
H c is much larger than possible energy losses. 

In this connection, one should mention a problem of energy balance in galactic halo 
gas, still far from being understood well. It is commonly believed that clustered SNe II in 
OB associations are the source which replenishes energy losses in the halo. Detailed calculations 
show however, that special conditions are needed to form vertical tunnels able to conduct energy 
from the disk to halo. It is shown in [39] that hot gas from a clustered SNe explosion can reach 
heights z ~ 1 — 3 kpc only if ~ 100 SNe exploded at distance z = 80 pc from the plane. Definitely, 
such events are rare in the Galaxy. Quite recently additional results have been obtained which 
exacerbates the problem. Based on recent all sky survey of HI gas, soft X-ray radiation, high 
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energy 7-ray emission and the 408 MHz survey, it has been found in [28] that radial distribution 
of halo gas has scale length Rh ~ 15 kpc. At the same time, galactic supernovae have different 
radial distributions: either exponential with radial scale length ~ 4 kpc for SNe I, or a ring 
with radius ~ 5 kpc and width ~ 4 kpc for SNe II. One might assume that SNe explosions in 
galactic central regions eject hot gas, which then fills the halo in radial direction over galactic 
scales R ~ 10 — 15 kpc. However, characteristic time for hot ejected gas to expand in radial 
direction out of solar circle is of order ~ 10 8 yr, which is comparable to its radiative cooling time 
3 — 10 x 10 7 yr. At the same time, the problem can have a simple solution if transfer of kinetic 
energy of evaporating material of DM-clouds to thermal energy is taken into account. Note, 
that this mechanism can be comparable (and out of solar circle can exceed) to SNe explosions, 
even if the total number of DM-clouds is smaller than the value accepted in previous sections 
- the rate H c decreased by factor of 10 is still of the same order as the SNe energy input 
rate. (This circumstance can be important if subsequent observations will restrict from above 
such parameters of clump uscules as M, M c , v c .) It worth to stress that since heating due to 
evaporating clumpuscules is much higher than radiation losses, it can support energy balance 
in halo on radial scales larger than vq: clumpuscules with radial scale rp ^ 10 kpc can maintain 
radial distribution of halo gas with R^ ~ 15 kpc (though r > 10 kpc seems more realistic). 

8. Discussion and conclusions. 

The effects connected with mass loss by DM-clouds due to their collisions with each other 
and photoevaporation by ionizing radiation, are interesting not only by possible observational 
manifestation, but also by their impact on dynamical and chemical evolution of the Galaxy. We 
discuss here only qualitative aspects of such influence - detailed analysis requires a separate 
consideration. As shown in §7, the mass loss rate by a single cloud ionized by the background 
UV radiation is M c ~ 10~ 13 M Q yr -1 . Thus, a clumpuscule of ~ 10~ 3 M Q loses all its mass 
on the Hubble time. Moreover, the life time of DM-clouds can be even less than the Hubble 
time since they cross regions near hot stars with enhanced UV flux. In addition, it is almost 
certain that the background UV radiation can be more intense in young Galaxy due to violent 
star formation. If one assumes that all evaporated material is trapped by inner parts of the 
Galaxy (R ^ 15 kpc), it follows inevitably that the total mass of luminous matter in the Galaxy 
must be comparable with (or higher than) that contained now in DM-clouds. Such a conclusion 
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would obviously contradict observations, if mean surface density in the Galaxy due to DM- 
clouds is order of £ ^ 10 2 M Q pc -2 as estimated in [18]. One can escape this paradox by 
limiting from above the cloud radius R c . Assuming clumpuscules to be gravitationally bound 
[6, 18], and with accounting Eq. (13) we obtain for cloud life time against photoevaporation 
t ev ~ 1.2 x W 23 T C R~ 1/2 . Thus, clouds survive in the Hubble time t H if R c < 1 - 3.5 x 10 12 
cm for cloud temperature T c = 2.7 and 5 K, respectively. (Note, that the restriction is weaker 
for higher temperature of DM-clouds: it scales as oc T 2 ). One should stress, however, that the 
evaporating gas (as well as dispersed in cloud collisions) can in principle settle down in the outer 
Galaxy and form an extended ionized disk, as suggested first in [40]. 

A lower limit for cloud radius R c can be obtained if possible interaction of DM-clouds 
with the Solar system is taken into account. The frequency of DM-clouds crossing the Solar 
system is 

v s ~ ttR 2 s Mv c , 
where M ~ //7i\R 2 r , Rs ~ 40 AU. This gives 

V vsro 

One can assume that v$ < (100 yr) _1 — more frequent collisions might result in detectable 
perturbations of planetary orbits. Then, for ro = 10 kpc and adopted parameters we get 
R c > 10 11 cm. 

The net mass per unit area deposited by evaporating clumpuscules can be found as 

A£ = £[(l + 2Wt CT ) 2 -l], 

where £ is the present surface density in clumpuscules. Requiring AS < £„, where T, v ~ 50 M Q 
pc -2 , we arrive at 



[(1 + 2t H /t ev ) 2 - 1] ' 

For the lower limit R c = 10 11 cm this gives £ < £„. This estimate conflicts with the lower 
limit, £ > 100 Mq pc -2 , obtained in [18]. The contradiction can be avoided if the covering 
factor / is two times smaller than estimated in [16] from the flux monitoring data. Note in this 
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connection that in [18] / ~ 10~ 4 is mentioned as a conservative value. [One should stress that 
since 2tn/tev < 1 arid t ev oc R c 1//2 the obtained upper limit weakly depends on the parameters: 
approximately as {fvs r o/ v c) l ^ A , though it is sensitive (inversely proportional) to the fraction of 
photoevaporated mass confined by the Galaxy] 

Transverse velocity of clouds should decrease proportionally v c ± oc R c to keep the 
duration of extreme scattering events constant. This means in turn, that orbits of DM-clouds 
are strongly stretched, so that the radial component (< 500 km s _1 ) is much bigger than the 
transverse one. 

With such restrictions the current mass rate processed in cloud collisions, vgM c oc 
f 2 R~ x , can reach about 2 — 20 M Q yr _1 and leads to overproduction of mass in the Galaxy. 
However, due to strong dependence on / it can be as small as < 0.4 — 4 M Q yr _1 for / < 10~ 3 . 

Although at present the origin of DM-clouds is far from being understood well, one 
may say tentatively that they were born on prestellar stages of the universe, and are thus the 
most old objects [13]. This means that they might have primordial chemical composition, and 
being mixed with the interstellar gas decrease its metallicity. In other words, models of chemical 
evolution of the Galaxy which incorporate input of mass from photoevaporating and collisionally 
dispersed DM-clouds, should predict less efficient enrichment by metals than those without DM- 
clouds. It follows, in particular, that since clouds concentrate to central regions of the Galaxy, it 
should make radial gradient of metallicity much softer than the observed one, or even to inverse 
it, (unless r is large enough, r > 10 — 15 kpc) - otherwise, star formation rate must increase 
strongly in central regions to balance mass input with primeval composition. 

We summarise our findings: 

1) Collisions of dense gaseous clumpuscules responsible for the extreme scattering of 
quasar radio emission, are frequent enough: depending on their total number in the Galaxy 
around 10 3 to 10 4 collisions per year can take place. 

2) After collisions gas of clumpuscules is heated and then dissipates and mixes with the 
interstellar gas. The corresponding mass input rate can reach 1 — 10 M Q yr _1 . In this case their 
impact on dynamical and chemical evolution of the Galaxy can be important. In particular, it 
might result in softening of radial gradients of the stellar metallicity. 
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3) On post-collisional stages shocked clumpuscules and their fragments appear as tran- 
sient optical sources with temperature and luminosity close to those of massive stars, the duration 
of these stages is about 10 days. 

4) On later stages expanding fragments can be seen in 21 cm, with duration of 1 Myr, 
and the total mass of such fragments in the Galaxy luminous in 21 cm is of 3 x (10 5 — 10 6 ) M . 
This gas can be distinguished from the interstellar HI gas by its high velocity dispersion. 

5) DM-clouds surrounding hot stars form HII halos around them with electron density 
of n e ~ 0.65 cm -3 , and can be seen in H a with emission measure EM ~ 20 cm -6 pc _1 even 
for stars far from the galactic plane where interstellar gas is diluted. 

6) Ionized gas outflowing the atmospheres of DM-clouds can be a substantial mass 
source - with the rate 2 — 20 M Q yr _1 - for the interstellar gas. If most of this evaporated gas 
is trapped by the Galaxy, clumpuscule radius can be constrained as R c < 3.5 x 10 12 cm from 
requirement the mass deposited to be less than the total luminous mass. In turn, it leads to the 
conclusion that orbits of clouds are strongly stretched with small transverse velocity component. 
In this case, the surface mass density in DM-clouds can be as small as 50 M Q pc -2 . 

7) In general, gas lost by DM-clouds has high velocity relative to the ISM (< 500 km 
s _1 [18]). Subsequent mixing with the interstellar gas transfers its kinetic energy to heat and is 
an additional heating source of the ISM. This can explain large radial scale length of the halo 
gas (~ 15 kpc) if the scale length of DM-clouds distribution (1) is high enough ( £ 10 kpc). 

Possible detection (or, contrary, non-detection at given threshold of sensitivity) of optical 
emission, including H a , and emission in 21 cm from clumpuscules will be of great importance for 
understanding the nature of dark matter and for more confident estimates (or firm constraints) 
of the contribution of dense molecular clumpuscules to the galactic mass. In principle, one can 
expect that dense atomic clouds, suggested in [11] as possible dark matter objects, also may 
manifest themselves in H a and 21 cm emissions through photoionization by UV background 
and collisions. However, their parameters (sizes, densities, velocity dispersion) are too uncertain 
now to make definite conclusions. In [7, 8, 12, 14, 41] baryonic dark matter is assumed to form 
aggregates similar to globular clusters - baryonic dark clusters - with a fraction of baryons in 
gas clouds (see for detailed discussion [12, 14]). Stability of such clusters implies that gas clouds 
are collisionless, however, even rare collisions might provide substantial mass supply into the 
ISM if the clouds are numerous enough. Also, photoevaporation due to external UV photons 
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must be important for such clouds, however, effects from neighbouring (brown dwarf) stars seem 
to be more influential [14]. 
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